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Abstract 

• We simulate two variants of quenched twisted mass QCD (tmQCD), with degenerate Wilson 
quarks of masses equal to or heavier than half the strange quark mass. We use Ward 
identities in order to measure the twist angles of the theory and thus check the quality of the 
tuning of mass parameters to a physics condition which stays constant as the lattice spacing 



is varied. Flavour symmetry breaking in tmQCD is studied in a framework of two fully 
twisted and two standard Wilson quark flavours, tuned to be degenerate in the continuum. 
Comparing pseudoscalar masses, obtained from connected quark diagrams made of tmQCD 
and/or standard Wilson quark propagators, we confirm that flavour symmetry breaking 
effects, which are at most 5%, decrease as we approach the continuum limit. We also 
compute the pseudoscalar decay constant in the continuum limit, with reduced systematics. 
As a consequence of improved tuning of the mass parameters at j3 = 6.1, we reanalyse 
our previous Bk results. Our main phenomenological findings are r^fx = 0.421(7) and 
B K = 0.735(71). 



1 Introduction 



Bk, the bag parameter of neutral K-meson oscillations, has been computed with sev- 
eral discretizations of lattice fermions. Until recently, the quenched Wilson fermion 
results of Bk were the least accurate, due to a limited control of those systematic 
sources of error, which arise form the lack of chiral symmetry in the regularization. 
This trend has been reversed in ref. [1], thanks to the implementation of twisted 
Wilson fermions [2] . The simulations of ref. [1] , besides introducing some novelties 
in the computation of Bk (twisted mass QCD (tmQCD) regularization, Schrodinger 
functional renormalization and RG running) have also been extensive from the com- 
putational point of view. In particular two tmQCD variants of the fermion action 
have been implemented (with twist angles ir/2 and 7r/4) at several inverse gauge 
couplings 0. Such a large collection of data enables us to address, in the present 
work, several other issues related to the tmQCD formalism, in the region of strange 
quarks. Clearly, the fact that we work in the quenched approximation is a limitation 
of the scope of the present work. 

In sect. [2]we present the details of our formalism. In our first tmQCD variant we 
introduce two twisted flavours, with twist angle n/2, and two standard (untwisted) 
flavours. In the second variant we only have two twisted flavours with twist angle 
7r/4. The tmQCD lattice action is 0(a) improved. All dimension-3 operators (cur- 
rents, scalar and pseudoscalar densities) are improved by introducing 0(a) Symanzik 
counterterms. This is essential to 0(a) improvement in the 7r/4 case and to those 
quantities of the tt/2 case which are not exclusively composed of fully twisted quarks. 
No improvement of the four-fermion operators is attempted. 

In sect. El we first of all examine the quality of the tuning of the mass parame- 
ters, fixed so as to ensure that the twist angle is equal to a target reference value (in 
our case tt/2 or vr/4). This, together with the requirement that all quark masses be 
degenerate and phenomenological quantities be measured at a reference pseudoscalar 
meson mass, constitute our constant physics requirement. It has to be maintained 
as we increase the theory's UV cutoff (i.e. approach the continuum limit). The twist 
angle is then measured with the aid of Ward identities and compared to its target 
reference value, used in the mass parameter tuning. We see that 0(a 2 ) cutoff effects 
are responsible for statistically significant discrepancies, which are nevertheless only 
a few percent. 

In sect. HI we measure the pseudoscalar masses and decay constants in both 
ir/2 and 7r/4 setups. In the former case we are also able to monitor flavour break- 
ing effects, arising from the twisted mass term in the lattice action. This we do by 
comparing the pseudoscalar masses and decay constants for mesons composed exclu- 
sively of twisted or untwisted valence quarks, as well as those made of one twisted 
and one untwisted quark. For the pseudoscalar masses, we find that statistically 
significant effects, which are nevertheless only a few percent at the coarsest lattices, 
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disappear as the lattice spacing is decreased. Thus flavour symmetry appears to be 
restored in the continuum limit. 

Our best results for the K-meson decay constant are based on a tmQCD Ward 
identity. They are free of the usual systematic uncertainties arising from current 
normalizations and improvement (i.e. ambiguities in the values of Z A , Z v and c A ). 
Extrapolating ir/2 and 7r/4 results to a common continuum limit gives a result in 
full agreement with earlier estimates. 

The more detailed analysis concerning the accuracy of the tuning of quark 
masses and twist angles, presented here, was performed after the publication of 
ref. [1]. The quality of the tuning was found to be satisfactory in all cases, save for 
the simulations at /3 = 6.1, mainly in the 7r/4 case. This is signalled e.g. by relatively 
large differences between the value of the target twist angle, set to ir/2 or 7r/4 through 
the tuning of the bare (subtracted) mass parameters, and the value obtained by 
computing the twist angle with the PCAC quark mass instead of the subtracted 
quark mass. The reason for this behaviour has been traced back to the value of K cr 
taken as input from the literature. Indeed, for an accurate determination of n CI it is 
crucial to fix its 0(a 2 ) ambiguities by following a constant physics condition in the 
approach to the continuum limit. Instead, the value of k ct ((3 = 6.1) quoted in [3] 
comes from an interpolation of data obtained from a constant physics condition at 
other values of (3. While the effect of relaxing the constant physics requirement was 
found to be negligible for the data of [3] , its impact on the tuning of twist angles is 
large. This is discussed in detail in Appendix lAl 

The (3 = 6.1 critical point has been hence determined afresh, and (3 = 6.1 
simulations with new mass parameters have been performed. Thus all results in the 
present work are generated from the datasets of ref. [1], except for those at (3 = 6.1, 
which are completely new. The new Bx((3 = 6.1) results can be found in sect. 
They induce a reanalysis of the continuum limit extrapolation of this quantity. Also 
in sect. we collect our detailed results of the kaon-to-pion four-fermion operator 
matrix elements, involved in the AI = 1/2 rule. Strictly speaking, these results are 
not physical, as they refer to four degenerate quarks, with masses close or above 
half the strange quark mass. They have been used, however, in ref. [4] in order to 
obtain the relevant four-fermion operator renormalization with Neuberger fermions, 
through a matching procedure of RGI matrix elements computed from both tmQCD 
and Neuberger regularizations. 

For the sake of legibility, all tables containing our results have been gathered in 
Appendix [El 
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2 General tmQCD formalism 



Twisted mass QCD has been designed to eliminate exceptional configurations in 
(partially) quenched lattice simulations with light Wilson quarks [2]. In its original 
formulation, it describes a mass-degenerate isospin doublet ij) of Wilson quarks for 
which, besides the standard mass term, a so-called twisted mass term ifiqipj^T 3 ^ is 
introduced. The properties of tmQCD have been studied in detail in [2], where, in 
particular, its equivalence to standard two-flavour QCD has been established^- We 
discuss here the main characteristics of this formulation, extended to more flavours, 
in ways analogous to those discussed in [8] and [1] . 

It is convenient to formalize our variant of tmQCD in terms of a twisted and an 
untwisted isospin doublet, denoted as ipt w = (V>i, $2) and = (^3, ^4) respectively. 
All flavours will eventually be tuned to be degenerate. The twisted (untwisted) 
isospin doublet is regularized in the standard tmQCD- Wilson (plain Wilson) fashion: 

S F = a A ^2 [ Aw(x)(D w + m ,tw + H5T 3 Mo,tw)'0twO»O 

X 

where D w is the standard Dirac- Wilson fermion matrix (with a Clover term) and r 3 
the Pauli isospin matrix. In this work, the Wilson plaquette action is the regular- 
ization of the pure gauge sector of the theory. For the rest of the notation, relating 
tmQCD to standard QCD (concerning field rotations, mass transformations etc.) 
see [1]. Here, what we are mostly interested in are the expressions for the renormal- 
ized quark mass in the twisted quark sector, given by the combination of standard 
and twisted mass parameters 



M R , tw = ^4 jtw + / 4 itw , (2.2) 
and the twist angle, defined in terms of renormalized masses as 

tana = ^ R,tw . (2-3) 

Also standard is the relation between renormalized and bare quark masses: the 
subtracted (unrenormalized) quark mass for Wilson fermions in denoted by om q = 
1/(2k) — 1/(2k ct ), k being the hopping parameter (2k = [omo + 4] _1 ). Whenever 
we need to identify the quark doublet / (with / = tw,w), we will denote the 
corresponding quantities by am q j and Kt . In a Symanzik 0(a) improved framework, 
the renormalized quenched quark masses for the untwisted flavours are given by 

m R _ w = Z m [m q>w (l + 6 m om q , w )] , (2.4) 



x for reviews on the subject see [5-7] 
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while the twisted quark masses renormalize as follows: 

WR,tw = Z m [m^ tw (l + b m am^ tw ) + b m anl tw ] , ( 2 - 5 ) 

MR,tw = ^Vo,tw(l + ^am q ,tw) • (2.6) 

In terms of the last two expressions, the twist angle may be expressed as [9] 

, / v a/x ,tw[l + (bp - b m )am qttw ] , . 

tan(a) = - — = . (2.7) 

Z P Z m [am^tw + b m (a/i 0: tw) 2 } 

We now discuss an alternative expression for the twist angle, obtained in terms 
of Ward identities in [9]. For the twisted sector of our theory, the PCAC and PCVC 
Ward identities read 

d„(A R )u,i2 = 2m R (P R ) 12 , (2.8) 
B v (V R ) u ,i 2 = 2ifi R {P R ) 12 , (2.9) 

and they are valid up to 0(a 2 ) for the Symanzik-improved operators [9] 

(P R )i 2 = Z P [l + 6 P am qitw ]P 12 , (2.10) 
{Ar)v,12 = Z A [1 + & A amq,tw] [A/,12 + ac A d u P 12 - ib A a^o M V v ^ 2 ] , (2.11) 
(Vr)„,12 = Z v [l + b v am qM ][V Utl2 + ac v d p T uPtl2 - ib v afi ,tw A u>12 ] , (2.12) 

where d v denotes the lattice symmetrized derivative and the operator subscripts 1, 2 
indicate quark flavours. By combining expressions (|2.8[) - (|2.12p with the standard 
definition of the PCAC bare quark mass [10] , 

d v \Av,vi + c A ad u P 12 ] = 2m tw P 12 , (2.13) 

the renormalized quark mass is computed aa^| 

Z A [1 + b A am qitw ] [ 2 r roi/l s 

Z P [1 + 6 P am q! twJ L J 

Finally, the PCAC expression for the twist angle, in terms of the above, is [9] 

, \ a/A),tw[l + (K + b P - 6 A )am q tw ] 

tan(a) = — : . (2.15) 

Z A [am tw + b A (afi QM ) 2 Z v } 

Yet another variant of the twist angle is obtained by expressing /i R tw 111 terms 
of the PCVC Ward identity. We define the PCVC bare twisted mass fromEl 

d u V„,i2 = 2i/i tw P 12 , (2.16) 



2 Eq. (|2.13p is to be understood in terms of correlation functions involving these operator inser- 
tions. In the present work Schrodinger functional correlation functions with pseudoscalar boundaries 
are implemented. The same considerations hold for Eq. (|2,16[) 

3 The tensor term proportional to c v in Eq. (12. 12[) vanishes upon differentiation of the vector 
current (Vr)„,i2. 



4 



and, using expressions (j2.8l) - (|2.12p . obtain for the renormalized twisted mass 

Z v [l + 6 v a"i q ,tw] 



MR.tw 



1 — am tw Z v b^ 



Mtw • (2.17) 



Z P [1 + 6 P am qi tw] 

The PCVC expression for the twist angle, in terms of the above, is 

Z v a// tw [l + (by -b A - Z v b v )am qyt w] /n 1 oN 

tan(a) = ~ — ■. — - — . (2-18) 

Z A [am tw + b A (a^ tw ) 2 Z v } 

In the present work, expression (|2.T|) is used for tuning the bare mass parame- 
ters [Iqm and to fix the theory to a specific twist angle. One also needs the 
value of k ct for the determination of am„i . This is known from previous Schrodinger 
functional computations, based on the Clover-improved theory with Wilson (un- 
twisted) quarks]^ Once the bare parameters are thus fixed to satisfy Eq. (|2.7p . one 
may use Eq. (|2.15|) and Eq. (|2.18|) in order to obtain independent estimates of the 
twist angle, which differ from the target value by 0(a 2 ) discretization effects. This 
provides a measure of the systematic uncertainties related to the tuning of the twist 
angle H 

Finally, we specify the twist angles, following the two cases of [1]. In the first 
case, known as fully twisted theory, the bare parameters of the twisted quark doublet 
are tuned so as to ensure that a = ir/2. This amounts to tuning Eq. (j2.5[) so that 
?fiR,tw = 0. The untwisted doublet is also tuned (through an opportune choice of 
k w ), so that tor. w = /ip^tw up to 0(a 2 ); cf. Eqs, (|2.4p and (|2.6p . In the second case 
we switch off the untwisted doublet, keeping just two twisted flavours; the twist 
angle is set to a = ir/4. This amounts to tuning Eqs. (|2.5p and (|2.6p so as to have 
ffiR.tw = MR,tw The detailed expressions used for these tunings may be read off from 
ref.'[l]. 

In order to make contact with the physical results we will be presenting, we now 
identify the four generic quark fields ipk = 1, - ■ ■ ,4) with physical (if degenerate) 
flavours. There is a different identification according to the problem in hand. When 
we discuss pseudoscalar masses, decay constants and Bk, in the tt/2 case we identify 
the twisted doublet with ?/>tw = (u,d) and the untwisted one with -0 W = (s, c). The 
same quantities in the 7r/4 case are addressed in terms of a single twisted doublet of 
a strange and a down quark; i.e. ^tw = (s,d). As our simulations are quenched and 
mass degenerate, this is sufficient to model two valence quarks. In this way we make 



The condition we implement for fixing the twist angle is neither the so-called pion mass deter- 
mination, nor the PC AC one of ref. [11]. 

5 The (re)normalization constants (Z A , Z P etc.) as well as the improvement coefficients (c A , b A 
etc.) used in this work are also taken from previous Schrodinger functional computations. Their 
values are all gathered in Appendix A of ref. [1]. Note that Eqs. (A. 10) and (A. 11) of that Appendix 
contain misprints; they should read b A = 0.086Cf<?o an( i ^v = 0.074Cf<?o- 
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full contact with the notation of the earlier tmQCD simulations of ref. [1]. For the 
results related to AS = 1 four-fermion operators, the physical flavour identification 
is more complicated. The reader is referred directly to ref. [4], where the issue has 
been addressed in detail. 



3 Cutoff effects of the twist angle 

We now turn to the determination of the twist angle from PCAC and PCVC rela- 
tions. The simulation parameters are gathered in Table [T] for the n/2 theory and 
in Table [2] for the 7r/4 one. The data are those of ref. [1], except at (3 = 6.1. As 
mentioned in the introduction, the run had to be repeated at this coupling, for rea- 
sons which will be discussed in detail below. The physical regime targeted in the 
runs of ref. [1] is that of the iT-meson, composed of two degenerate valence quarks. 
As detailed in that work, in practice this means that the choice of quark masses is 
such that the K- meson in the tt/2 theory is computed in the range 640-830 MeV 
and extrapolated to the physical point at 495 MeV (i.e. roM]^ hys = 1.2544). In the 
7r/4 theory we are instead able to simulate with quarks corresponding to a physical 
kaon of about 495 MeV. The only exception is the (3 = 6.45 case, in which extrap- 
olations from higher mass values were the only option, as simulations with quarks 
corresponding to a physical kaon require prohibitively large lattice sizes. As the 
present work is based on the same runs, the same physical point for all four flavours 
is targeted here. 

Following the discussion of sect. [2j the bare mass parameters (i.e. standard 
hopping parameter(s) and twisted mass /io.tw) are tuned at each (3 value so as to 
keep the quarks degenerate and the twist angle fixed at a = 7r/2,7r/4. The n cr 
estimates used in the bare parameter calibration are taken over from the literature; 
see Table EJ Those of ref. [12], as well as the one provided to us by the ZeRo 
Collaboration^!, are the result of a direct computation of the PCAC Ward identity at 
the corresponding bare coupling, using 0(a) improved (untwisted) Wilson fermions, 
whereas the estimates from ref. [3] are the result of an interpolation in (3 of the data 
of ref. [12]. We will see that this k ct estimate at (3 = 6.1 results in poor tuning of 
the bare tmQCD parameters. Thus, following the procedure described in ref. [12], 
we have recomputed k ct at (3 = 6.1, finding a different value. This has been done 
by an independent run, performed at lattice volume L/a = 16, T/L = 2 and four 
PCAC quark mass values in a range 0.025 < am av < 0.05, similar to that of ref. [12] 
(where m av is defined). Our k ct result is obtained by linear exrapolation in the 
PCAC quark masses, ensuring that our systematics resemble those of ref. [12]. 

All observables computed in this work are obtained from the large time asymp- 

6 We thank I. Wetzorke for providing us with the K CT (f3 = 6.3) value, obtained in the context 
of ref [13]. 
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totic limit of operator correlation functions with Schrodinger functional (SF) bound- 
ary conditions. The notation is standard, following closely that adopted in e.g. [1]. 
For instance, /a, 12 denotes the Schrodinger functional correlation with a fermionic 
operator Aq^ in the bulk and a time-boundary pseudoscalar operator C275C1 at 
x = 0. All such correlation functions are properly (anti)symmetrized in time, when 
used to extract quark masses, effective pseudoscalar masses and decay constants. 
The bare PCAC and PCVC quark masses of Eqs. (|2.13|) and (|2.16p are obtained 
from the ratios 



From these and Eqs. (|2.14p . (|2.17p we compute the renormalized masses oma,tw 
and a/iR,tw! from Eqs. (I2.15p . (|2.18|) we compute cot (a). These results are reported 
in Table [H The errors are statistical; we have checked that systematic errors due 
to the uncertainties of K cr , the (re) normalization parameters and the improvement 
coefficients are an order of magnitude smaller than statistical ones. 

Let us comment on our ir/2 results first. The two a/XR tw values (computed from 
Eq. (12. 6p and Eq. (j2.17D ) have a small but statistically significant discrepancy of less 
than 3%. Similarly, the result for amRtw is also statistically different from the target 
value aniR^ = 0. The two estimates of the twist angle are only a few percent off 
the target value n/2, but the small discrepancies in the a//R >tw evaluations do not 
carry over to the two results for cot (a), which are nearly always compatible within 
errors. All observed deviations from the expected target values are attributed to 
0(a 2 ) discretization effects. 

The 7r/4 results display largely the same characteristics, but the mass range 
spanned by the data (for (5 < 6.45) is too narrow to discern the details of the 
dependence of the twist angle upon the quark mass parameters. What we see is 
that the difference between the two twisted mass estimates /iR,tw an d the standard 
mass estimate m^ tw is a small but statistically significant effect. The same is true 
of the two Ward identity estimates of the twist angle, which also differ from the 
target value 7r/4 by a small amount. 

Finally, we study to which extent the Ward identity estimates of the twist angles 
approach their target values 7r/2,7r/4 in the continuum limit. We do this by first 
computing, at fixed value, the quantity cot (a) at the kaon mass reference scale 
Mk = Mj^ hys = 495 MeV. This is done by extra/interpolating our data as a function 
of the pseudoscalar effective mass-squared, expressed in physical unitslll The result 
thus obtained is plotted against (a/ro) 2 in Fig. [TJ A linear extrapolation of cot (a) 

7 The pseudoscalar mass in question is obtained from a correlation function consisting exclusively 
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to the continuum limit turns out to be unreliable (very large x 2 /d.o.f.), as the data 
do not display a monotonic behaviour, with fluctuations which are much larger than 
their errors. In any case, these fluctuations are only a small effect, reflecting the 
overall uncertainty of the tuning of the twist angle to a constant target value (which 
amounts to a condition of constant physics as we approach the continuum limit). 
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Figure 1: Lattice spacing dependence of cot (a), computed at the physical kaon mass 
point, for twist angles ir/2 (circles) and 7r/4 (squares). The errors are smaller than 
the points. 



4 Flavour symmetry breaking, meson masses and decay constants 

In order to investigate flavour breaking effects, we compare effective pseudoscalar 
masses, obtained from correlation functions with different combinations of twisted 
and standard Wilson quark propagators, on our quenched configuration ensemble of 

of twisted quark propagators. In the next section, this effective mass is denoted as (roM^f ) for the 
7r/2 case, while for the 7r/4 one we use (roiWjJ). 
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tmQCD at a = tt/2. The rationale behind these computations is as follows: we are 
dealing with a (quenched) lattice QCD model with four degenerate quark flavours, 
two of which are twisted. The massless continuum theory has an SU (4)l <8> SU (4)r 
chiral symmetry with 15 degenerate pseudoscalar Goldstone bosons; with massive 
fermions we are left with the vector flavour symmetry SU(4)y. In the regularized 
theory, the Wilson term breaks the symmetry induced by axial transformations, even 
in the absence of quark mass parameters. In twisted mass QCD some of these axial 
symmetries are interpreted as part of the SU(4) flavour symmetry. With our choice 
of twisted mass terms one then finds that on the lattice the SU(4) flavour symmetry 
is reduced to the subgroup U(l) 0(7(1)0 SU(2). Thus, at finite lattice spacing, we 
expect that the "twisted" charged Goldstone bosons, the "untwisted" ones and the 
"twisted- untwisted" ones will differ in mass by terms which are like 0(a 2 ) (recall 
that we work with a Symanzik-improved action). Based on this approach, we have 
provided a first summary of our findings on this symmetry breaking in [1]; here 
we present our full results. The same approach for monitoring such flavour break- 
ing effects has also been implemented (with an action with two 7r/2-twisted isospin 
doublets and no Clover term) in [14]. Both works focus on pseudoscalar masses in 
the kaon region. For similar results closer to the chiral limit (for a single (3 value), 
see [15]. It is fairly straightforward to measure these flavour breaking effects, as the 
corresponding correlation functions involve only connected diagrams. Recall that 
there is also a flavour breaking effect between the charged and neutral "twisted" 
pseudoscalar, which is harder to monitor, as the measurement of the neutral pion 
mass involves disconnected diagrams; see [16] Studying this flavour breaking is be- 
yond the scope of the present work. 

We have measured pseudoscalar effective masses aM??(xo), using suitable time- 
dependent correlation functions /( Ar ). j (xo) (with i,j distinct flavour indices) of the 
temporal component of the renormalized axial current (Aji)ij 



/(A R ) tJ fro ~ a) 

/(A R )y(z +a) 



(4.1) 



These quantities are suitably (anti) symmetrized in time and averaged over plateaux, 
as detailed in [1]. In the language of the twisted action Eq. (|2.ip at twist angle vr/2, 
the corresponding lattice correlation functions to be used in Eq. (|4.ip are: 

* f(An)sA x o) -> ^[/(A R ) a >0) -«/(v R ) sd (zo)] , 

* f(A R ) ud ( x o) -> -«7(v R ) ud (^o) , 

* /(A R )«(»o) -» f(A R ) sc ( x o) ■ 
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The first correlation function is composed of a tmQCD Wilson quark propagator and 
a standard QCD Wilson quark propagator. Prom it we derive the "-fT-meson" effec- 
tive mass, denoted by aM^[ . The second correlation function, composed exclusively 
of the tmQCD Wilson quark propagator, provides the charged "pion" effective mass 
aM^j. Finally, the third correlation function, composed exclusively of the standard 
QCD Wilson quark propagator, provides the "D s -meson" effective mass aiW™. Since 
all quark masses are tuned to be degenerate, these are three of the 15 degenerate 
Goldstone bosons in the continuum limit of our quenched theory. At finite lattice 
spacing tmQCD induces flavour breaking discretization effects, which are monitored 
by comparing the values of the three effective masses. 

The corresponding currents, inserted in the above correlation functions, are the 
following Symanzik-improved quantities, taken over from Appendix B of [1]: 

(A-a.) V7 ad = Z A [1 + -b A amq >w ][A UjSd + ac A d u P sd - i-b A a/j, 0M V u>s d] , (4.2) 

(VR)u,sd = Z v [l + -Kam^ w ][V v>sd - i-b v a/j, 0j twA u>s d] , (4.3) 
{VR) UyU d = Z Y \y v>ud - ibya^twA^ud] , (4.4) 
(A R ) VtSC = Z A [1 + b A am qtW }[A UtSC + ac A d u P sc ] . (4.5) 

Again the 0{a) tensor-like counterterm in the above vector currents has been omit- 
ted, since it drops out in the correlation functions (it is a sum over space, with 
periodic boundary conditions, of a discrete spatial divergence). 

Another estimate of the effective pseudoscalar meson mass, denoted as aM^, 
is obtained by using the correlation f( Pn ) ud in Eq. (|4.ip . with 

(PR)ud = Z P [l + b P am qM ]P ud . (4.6) 



The results are collected in Table [5l The agreement between Mf^ and M^ d 

is excellent. In Fig. [2] the ratio (M sd /M sc ) 2 , plotted against (M sd /Mjf ys ) 2 with 

M^ ys = 495 MeV, is compared to (M ud /M sc ) 2 . Note that the statistical errors 
of these ratios turn out to be very small, due to the strong correlations between 
numerator and denominator. We see that at f3 = 6.0 the two ratios are incompatible; 
their deviation from unity (which quantifies flavour breaking) is at most a 5% effect. 
As we approach the continuum limit at (5 = 6.3, the two ratios become compatible 
and their deviation from unity reduces to 1 — 2%. Our conclusion is that flavour 
symmetry breaking effects at the mass ranges we are considering appear to be under 
control, diminishing fast as the continuum limit is approached. Besides this general 
conclusion, there are a couple of observations to be made: (i) At (3 = 6.0 we find 
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that M ud < M sd , in agreement to the findings of [14] (see Fig. 13 of that work)H; (ii) 
at fixed reference mass (M sc /M^- hys ) 2 , these mass ratios do not display a monotonic 
dependence on j3. This is probably a small cumulative effect of the many systematic 
uncertainties of the mass tuning procedure. 

A similar analysis is performed for the pseudoscalar meson decay constants. 
In the Schrodinger functional framework they are obtained from the axial current 
correlation functions / Ar , properly normalized by the boundary-to-boundary corre- 
lation function /i; see [17] for details. For the tmQCD tt/2 case under investigation, 
the specific expressions are (in the large time asymptotic regime): 

F sd « V2{Mf d L 3 )-^exp[(x -T/2)Mf ] AarW^o) - »/(v R w(^o) (4J) 

V Jl,sd 

F ud « 2{Mf d I?yy* exp[(x - T/2) M f d ] , (4.8) 

V Jl,ud 

F sc « 2{M«I?)-W exp[(x - T/2)M 5 f ] , (4.9) 

V ■' 1 > sc 

The quantities F ud and F sc are obtained from these expressions in a range of 
xq 3> in which the pseudoscalar effective masses have been extracted. 

A second method for computing F ud is based on the PCVC relation Eq. (|2.9p . 
expressed in terms of Schrodinger functional correlation functions: 

- M ud f iVR)ud (xo) = 2in RM f {PR)ud (x ) . (4.10) 

The corresponding decay constant is computed as 

^ w - 4 ff^ L3 )" 1/2 ex PK*° - T / 2 )<"] /(PR ^ 0) ■ (4.11) 

The results for the decay constants are collected in Table EJ Note the excellent 
agreement between the F ud and F ud results, at all j3 values, which is analogous to 
the agreement between the two cot (a) computations presented in sect. O In Fig. [3] 
the ratio F sd /F sc is compared to F ud /F sc . The situation is qualitatively analogous 
to that of the mass ratios presented above, in that flavour breaking effects tend to 
vanish as the continuum limit is approached. These effects range form 13% to 3% 
with increasing f3; however we show below that these estimates depend heavily on 
the improvement coefflcent of the axial current. 

The above results have been obtained with the normalization constants and 
improvement coefficients determined in various ALPHA Collaboration publications. 



The authors of [14] call M K + what we call M u d and "kaon with Wilson strange" what we call 
M, d . 
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In order to monitor their influence on our data, we have repeated the analysis using 
the LANL Collaboration results for the same quantities (for details, numerical values 
and references see Appendix A of [1]). The comparison of these results, extrapolated 
to the physical kaon point, is displayed in Fig. [H From it we draw the following 
conclusions: 

• The ALPHA and LANL results for each of the decay constants are compatible, 
the only exceptions being F sc (incompatibility) and F s< i (near incompatibility) 
at (3 = 6.0. More detailed tests indicate that the main source of incompatibility 
lies in c A {j3 = 6.0). 

• The ALPHA and LANL results for F u d are fully compatible to those for F u d- 
All three are independent from the axial current (cf. Eq. (j4.8f) and Eq. (|4.1ip ) 
and thus from c A . Moreover, F u d does not depend on any normalization con- 
stants and/or improvement coefficients (cf. Eq. (|4.1ip ). 

• Unlike the ALPHA results, at (3 = 6.0 the three LANL estimates F sc , F s d and 
F U( i are fully compatible □ and show less scaling violations over the range of 
simulated couplings. This implies that the 0(a 2 ) discrepancy shown in Fig. [3] 
is significantly reduced if c^ ANL is used instead of c ALPHA . Thus the tmQCD 
flavour breaking effects under scrutiny, when monitored by decay constant 
ratios, appear to be obscured by the uncertainty in the c A determination. The 
mass ratios of Fig. are a more reliable monitor of tmQCD flavour breaking. 

The continuum limit estimates for the various decay constants of the ir/2 
case are obtained by linear extrapolation in (a/ro) 2 . Strictly speaking, Symanzik- 
improved quantities such as the decay constants, contain some improvement coef- 
ficients which are only known in perturbation theory (cf. b A and by). This means 
that there are also 0(o<7q) discretization errors. We have explicitly checked that the 
influence of the corresponding counterterms is negligible in practice and therefore 
the dominant discretization error is indeed 0(a 2 ). 

From Table [6] we see that our continuum limit results are compatible across 
all flavour combinations considered. Moreover we note that they do not change 
substantially if the point of the coarsest lattice ((3 = 6.0) is removed. From Fig. 0] 
we see that the continuum limit extrapolations with ALPHA and LANL data are 
also compatible. We have also confirmed that this conclusion remains valid if the 
(3 = 6.0 data are included in these extrapolations. 

9 The (3 = 6.0 LANL values at the physical kaon mass are r F sc = 0.4205(63), r F sd = 0.4275(74) 
and r F ud = 0.4292(67). 

10 This is not in accordance with the findings of ref. [18] for F aCl the continuum limit of which was 
obtained without the /3 = 6.0 result. This difference is explained by the bigger statistical sample 
(0(1000) configurations) and the simulations down to a finer lattice spacing (/3 = 6.45) of that 
work. 
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Figure 4: Pseudo-scalar meson decay constants, obtained with ALPHA and LANL 
normalization and improvement coefficients at the physical kaon mass. The twist 
angle is ir/2. The results for ro-F n d (based on a Ward identity) are also shown. The 
continuum extrapolation (linear in (a/ro) 2 ) is obtained without the (3 = 6.0 data. 
The ALPHA and Ward identity abscissae are slightly displaced for clarity. 
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We now pass to the computation of the kaon masses and decay constants in the 
7r/4 case. The corresponding lattice correlation function to be used in Eq. (|4.ip is: 

• /(A R )„,(zo) ^|[/(Aa) sd (^o)-i/(v R ) sd (^o)] , 

where now flavours s, are both twisted. From it, the estimate Mf? is obtained. A 
second estimate Mf5 can also be obtained from the pseudoscalar density correlation 
function /( PR ) ad , with 

(P R ) sd = Z P [1 + & P am q ,tw]P s d • (4.12) 
The decay constant is obtained from the correlation function 

F sd « v^Mjf L 3 )-V2 eX p[ (;Eo _ T/2)M s c f] /(VR)sd(X0) ~!j VR)sd(X0) - . (4.13) 

An alternative derivation is based on the continuum PCAC relation 

m (^r)*«j = 2M Ritw (P R ) sd = 2v / 2MR,tw(i :, R) s d • (4.14) 

The last equation is derived by taking into consideration that Eq. (|2.2p reduces to 
Zp y/2^o,tw for degenerate quark masses wr^w = MR,tw We thus obtain the decay 
constant estimate 

F sd « -4v / 2^(Mf L3)-V2 exp[(xo _ r/2)Mf] /(p ^ 0) . (4.15) 
M sd vh,sd 

The results for the decay constants, obtained with ALPHA estimates for the 
normalization constants and improvement coefficients, are collected in Table [7J In 
most cases there is again full compatibility between the masses vqM^ and tqM^ , 
as well as the decay constants vqF U( i and tqF uc i, at all (3 values. The comparison 
with the LANL results, made at the physical kaon mass, is displayed in Fig. [5] From 
it we draw the following conclusions: 

• The ALPHA and LANL results are compatible, except at f3 = 6.0. 

• Beyond (3 = 6.1, ALPHA and LANL results for roF s d are compatible to those 
for r$F S( i, the latter being independent of any normalization constants and/or 
improvement coefficients (cf. Eq. (14. 151) ) . As the ALPHA r^Fgd estimate scales 
like vqF sc i for the whole (3 range, the two have almost identical continuum 
limits. 

• Compared to the ALPHA and Ward identity results, the LANL ones display 
a better scaling behaviour over the whole range of simulated couplings. 
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The continuum limit extrapolations for the decay constants, linear in (a/ro) 2 , are 
also displayed in Table [7J We see that they are all in good agreement and do 
not change substantially if the point of the coarsest lattice {(5 = 6.0) is removed. 
The same extrapolations (for all [3 values) performed with the LANL data yield 
r oFsd = 0.404(6), which is incompatible to the ALPHA and Ward identity estimates, 
due to its small error. If however the j3 = 6.0 point is removed, we obtain the 
compatible result vqF S( i = 0.411(9), which is the situation displayed in Fig. [SJ 

The overall conclusion is that the continuum results are remarkably stable for 
the different flavour combinations, and different tmQCD regularizations (i.e. ir/2 
and 7r/4 cases). The best result for Fr in the continuum limit is obtained with a 
constrained fit of F u d (for the ir/2 case) and F s d (for the 7r/4 case). This choice is 
dictated by the absence of (re)normalization and improvement coefficients in these 
quantities, which amounts to the elimination of one source of systematic errors. Our 
final Fk estimate is 

r F K = 0.421 ± 0.007 (4.16) 

which agrees nicely with the previous ALPHA result r$Fx = 0.415±0.009 of ref. [18], 
and the ^LF-Collaboration one r F K = 0.410 ± 0.011 of ref. [19]. 



5 Four-fermion operators 

We now pass to the discussion of our tmQCD results concerning the matrix elements 
of the parity-odd four-fermion operators 

Qf = bPllnfa] [^37/^75^4] + [^17/^75^2] [^37/^4] ± [2 4] (5.1) 

between pseudoscalar states. After attributing appropriate physical quark flavours 
to the four-quark fields ipk(k = 1 5 ■ • • , 4), we compute the ratios 

(7T + \Qf\K+) 

± " {7r+\A OtUd \0){0\A o , au \K+) ' (5 ' 2) 

These ratios are the core quantities for the calculation of the kaon-to-pion weak 
matrix elements related to the AI = 1/2 rule. Note that our computations do not 
directly provide the physical K — > tt matrix elements of interest, as our pseudoscalar 
mesons (pions and kaons) are degenerate and at best as light as the physical kaon. 
Nevertheless our results have been used in ref. [4] in order to obtain the renor- 
malization constants of Qf with Neuberger fermions. This is achieved through a 
matching procedure involving the corresponding RGI operators, computed with Wil- 
son fermions. For details on the method and the notation, see ref. [4]. A different 
relabelling of the operator quarks ipk(k = 1 } • • • >4), allows the identification of the 
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Figure 5: Pseudoscalar meson decay constants, obtained with ALPHA and LANL 
normalization and improvement coefficients at the physical kaon mass. The twist 
angle is 7r/4. The results for TqF u ^ (based on a Ward identity) are also shown. The 
continuum extrapolation (linear in (a/ro) 2 ) is obtained without the /3 = 6.0 data. 
The ALPHA and Ward identity ordinates are slightly displaced for clarity. 



matrix element of Qf between pseudoscalar states with (K°\Q^ S 2 \K°). More- 
over, in the quenched approximation we have the identification Bk = (3/4)72.+ ; see 
refs. [1,20,21] for details. 

Our aim here is twofold: First, we wish to report our B k results at (3 = 6.1, 
computed with the new value of k ct . Second, we also list the values of 72 _ at all 
couplings; the analysis of ref. [4] was based on these results. As the continuum 
extrapolation of 72+ has also been performed in ref. [4], it needs not to be repeated 
here. Nevertheless, we will discuss in some detail the continuum extrapolation of 
Bk (which is that of 72+, since Bk = (3/4)72+), as this amounts to an update of 
our older Bk result of ref. [1]. This new Bk analysis has also been presented in 
ref. [6]. 

All computational details are identical to those of ref. [1]. Following that work, 
we construct three versions of the ratio 72+ , differing in the structure of the 0(a) 
bilinear counterterm operators in its denominator. The results for 72+ at (3 = 6.1, 
are collected in Table [HJ the corresponding RGI bag parameter Bk is shown, for all 
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Figure 6: Uncertainties on Bk related to the 0(a) improvement of bilinears. The 
left panel displays ir/2 data, while the right panel shows 7r/4 data. 



lattice spacings, in Tableland Fig. [6j The determination of Bk in the continuum 
limit involves a linear extrapolation in a, with the two different regularisations (vr/2 
and 7r/4) combined in a fit constrained to a common value at zero lattice spacing. 
It turned out that one of the most relevant sources of cutoff effects is related to the 
arbitrariness of the denominator 0(a) counterterms mentioned above. For instance, 
using either the values for Z A ,Z v ,c A determined by the ALPHA Collaboration or 
those obtained by the LANL group [22] results in sizeable effects on Bk at (3 = 6.0. 
At (5 = 6.1 we also discern discretization effects in the 7r/4 case (see Figured]). This 
signals the presence of large 0(a 2 ) ambiguities in Bx far from the continuum limit. 
Combined linear+quadratic extrapolation of the data proved to be unreliable, since 
the curvature of the quadratic term dominates the result also close to the continuum 
limit. Linear fits, excluding the (3 = 6.0 data, give 

B K = 0.668(45) , with c^ LPHA ; (5.3) 
B K = 0.737(51) , with c\ ANL ; (5.4) 
B K = 0.744(42) , with c A = ; (5.5) 



19 



whereas once also (3 = 6.1 is excluded we obtain 

B K = 0.735(71) , with c ALPHA ; (5.6) 
B K = 0.771(80) , with c\ mh ; (5.7) 
B K = 0.780(67) , with c A = . (5.8) 

In the ir/2 case, the error decrease in the extrapolations upon including the (3 = 6.1 
point (for which LANL and ALPHA data are fully compatible), is marginal. For 
instance, combined linear extrapolations with c ALPHA , excluding only the (3 = 6.0 
data for vr/2 and those at /3 = 6.0,6.1 for vr/4, yield B K = 0.739(66). 

The above results indicate that the extrapolation of the LANL data is the most 



stable. In spite of this, the ALPHA data is considered to be the best estimate, on 
grounds related to the systematic uncertainties in the derivation of Z A , Z w and c A , 
as explained at length in ref. [1]. Since the difference between ALPHA and LANL 
results is only significant at (3 = 6.0 and (3 = 6.1, we have conservatively discarded 
these data points in the continuum extrapolation, illustrated in the left panel of 
Figure [71 The final results are: 

I^r = 0.735(71) , (5.9) 
B^(2 GeV) = 0.534(52) . (5.10) 

When comparing with the result (Bk = 0.789(46)) quoted in [1], we must take 
into account that the revision of the [3 = 6.1 data has important consequences for 
the continuum limit extrapolation. In the analysis of ref. [1] (cf. right panel of 
Figure [7]), good scaling behaviour appeared to set-in rather abruptly at (3 = 6.1, 
with continuum limit extrapolation becoming stable only once the (3 = 6.0 points 
were discarded. On the contrary, the new (3 = 6.1 data interpolate in a smoother 
way those at [3 = 6.0 and (3 = 6.2 (cf. left panel of Figure [7|). This however, implies 
a worsening of the scaling behaviour, with continuum limit extrapolation becoming 
stable only upon discarding our (3 = 6.0 and /3 = 6.1 results, as detailed above. 

The value of Bk is shown in the bottom panel of Figure [7] alongside other repre- 
sentative results in quenched QCD found in the literature [23-30]. Our result is the 
only quenched result which has simultaneously eliminated the systematic uncertain- 
ties related to renormalisation (both at a reference scale and from the point of view 
of RG running), ultraviolet cutoff dependences, and finite volume effects (within the 
available accuracy). On the other hand, the control of the mass dependence of Bk 
with Wilson fermions is still not as accurate as with e.g. Neuberger or domain wall 
fermions. 

The 1Z- results for all [3 values are displayed in Tables [10] and [TT1 The extrap- 
olation of this data to the continuum limit has been presented in ref. [4]. 



20 



<x = jc/2 h 


— B — i 


a = ji/4 ^ 


-e— < 


cont. extr. t 


— A— i 



0.9 



<m 0.8 



0.7 - 



0.6 



T — 1 — | — 1 — 1 — 1 — 1 — | — 1 — 1 — 1 — 1 — | — 1 — 1 — 1 — 1 — | — 1 — 1 — 1 — 1 — | — 1 — 1 — 1 — I - 


1 1 


- K cr ((3=6.1) from this work 














1 
















- / 


i.' * " 












i i 


, , , , 


, , , , 


, , , , 


, , , , 







0.9 



0.02 0.04 0.06 0.08 0.1 

a/(2r ) 



0.7 - 



0.6 



T — 1 — | — 1 — 1 — 1 — 1 — | — 1 — 1 — 1 — 1 — | — 1 — 1 — 1 — I - 


III! 


11,1 

i 




- K cr (P=6.1) from [3] 


















M 






. "i 


i "-""="=" 












. . 


, , , , 


, , , , 






, , , , 




0.02 0.04 0.06 0.08 0.1 

a/(2r ) 



<m 0.8 




MILC03[21] Babich et al 06 [23] SPQ cd R 04 [24] JLQCD 97 [26] 
CP-PACS 01 [20] BosMar 03 [22] ALPHA 06 Lee et al 04 [25] 



Figure 7: Left: Continuum limit extrapolation of Bk- The tt/2 ordinates have been 
slightly displaced for clarity. Right: Same as before, but with Bk{(3 = 6.1) from [1]. 
Bottom: Comparison with other quenched results; different fermion discretizations 
are as follows (from left to right): [23,24] domain wall; [25-27] overlap; [28] Wilson; 
[29,30] staggered. The ALPHA point is that of the present work (Wilson tmQCD). 



6 Conclusions 

In this work we have completed our study of basic kaon weak matrix elements in 
quenched Wilson tmQCD. Our final value for Fk is the best controlled quenched 
result obtained with Wilson fermions from the point of view of systematic uncer- 
tainties. We have also provided a final value for Bk, with an error that, in our view, 
reflects faithfully the best accuracy that can be expected for this quantity in the 
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absence of full 0(a) improvement. Finally, we have performed a thorough study of 
flavour breaking in our version of tmQCD, confirming that it does not introduce any 
uncontrolled systematics in our results. 

The dominant source of uncertainty left in the quenched approximation (cer- 
tainly so for Bk) is related to the lack of full 0(a) improvement, which amplifies 
the error of the continuum limit extrapolation. Thus, if Wilson fermions are to be 
used in the future in the determination of weak matrix elements, the use of tmQCD 
variants that embody automatic 0(a) improvement [31] may prove essential. Two 
important aspects of the tmQCD approach are crucial in the context of weak matrix 
elements: the tuning of mass parameters, in particular of the twist angle, has to be 
controlled to high precision; and flavour symmetry breaking effects should be rea- 
sonably small, as in the present study. In conclusion, the present work demonstrates 
that, once the tuning of the twisted angle and flavour symmetry breaking are under 
control, tmQCD may offer a convenient alternative to other discretizations. As we 
are entering the era of tmQCD simulations with dynamical quarks [32], it will be 
important to explore these issues in future unquenched studies. 
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Appendix A The effect of an offset in the k ct estimate 

In this appendix we present the results for cot (a), the pseudoscalar effective mass 
and the decay constant, computed with quark mass parameters tuned with the 
value of n cr ((3 = 6.1) quoted in ref. [3]. As shown in Table O this value, obtained 
by interpolation over a range of re cr (/3), is roughly only 0.1% off the estimate of 
the present work, computed directly at (3 = 6.1. This apparently small offset is 
nevertheless a discrepancy of 15 standard deviations and has significant consequences 
in the quantities of interest, which we now examine in some detail. 

We denote by K cr the value computed on the constant physics condition of 
ref. [12]; at (5 = 6.1 this computation has been performed in the present work. The 
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estimate at (3 = 6.1, obtained through interpolation of data computed at several 
other gauge couplings [3], is parametrized ctS — ^cr 

[1 + 6 CT ] . Let us keep track 

of the effect of this offset in the tuning of the various quark masses, starting with 
the 7r/2 case. The tuning of the hopping parameter k w of the untwisted doublet is 
done in our simulations by requiring that the pseudoscalar meson made of the two 
untwisted flavours, has a mass M^? equal to a given value, fixed between 640 MeV 
and 830 MeV. This procedure is unaffected by the computation of purely 

untwisted quantities such as roF sc . What is affected by the offset 8 CT is our estimate 
of the subtracted masses, which now become 

am' qtf = — - -^j- = am^f + (6am CI ) , (/ = tw,w), (A.l) 

with 

(5am CT ) = S ct /(2k ct ) H , (A.2) 

where terms of 0(5 2 r ) have been dropped. So now Eq. (|2.4|) becomes 

am R, w = am R,w + Z m (5am cr ) H . (A.3) 

where O(a) counterterms are dropped, when multiplied by 5am cr . The requirement 
of quark mass degeneracy now reads m' Rw = ^Rtw' ^- e - ^ ne on?se t filters through to 
the twisted mass, which now becomes 

a /4,tw = aMR,tw + Z m (8am cx ) H . (A.4) 

This clearly induces an offset Zp 1 Z m (5am cr ) in the tuning of the bare twisted mass 
parameter, which we therefore denote as an' otw . Finally, the hopping parameter of 
the twisted doublet is tuned to a value Rt w , corresponding to 

om q ,tw = T^r = am q ,tw + (<5am cr ) , (A.5) 

by requiring the vanishing of the renormalized light quark mass, which is now written 

as 

am' RM = Z m [am qtw (l + b m afh' qM ) + b m {afi' 0tw ) 2 } . (A.6) 

It is important to note that the above quantity is not the true standard quark 
mass in the twisted sector. Considering that k ct (and not k' ct ) is by definition the 
reliably estimated critical point, the true renormalized quark mass, for the hopping 
parameter R tw , is expressed in terms of m q tw (cf. Eq. (|A.5p ) as am qitw = am^ — 
(5am cr ). It turns out to be non-zero: 

awR.tw = Z m [am q>tw (l + 6 m am qitw ) + b m (a^ 0M ) 2 } = -Z m (5am CT ) . (A.7) 
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The second expression has been derived by implementing the vanishing of Eq. ()A,6j) . 

The bottom line is that we now have a theory characterized by the bare param- 
eters k w , «twj /4w ( or ' equivalently, m qjW , m q]tw , /j,' tw ) which correspond to the same 
heavy quark mass m RjW , a small but non-zero light quark mass m Rjtw and a twisted 
mass MRtw ^ twist angle defined through the ratio of mass parameters a/4uw and 
am Rtw , is tautologically equal to ir/2. Instead, the true twist angle of the theory is 
given by 



, /-a °^R,tw am R,w am q , w 
tanfa ) = _ — = — — — = — - — - — I , (A. 8) 

which may differ significantly from the target value n/2. 

In Table \12\ we show the values of cot (a). It is clear that they are completely 
incompatible to the target value of ir/2 (and 7r/4; see below), on which the mass 
parameter tuning is based. This is in contrast to the small deviation from the target 
values of the cot (a) estimates given in Table HI which simply reflects the presence 
of 0{a 2 ) effects. 

In Table [T3l the results for the various pseudoscalar effective masses and decay 
constants are presented for the it/ 2 case. The quantities r^M^ and ro-F sc are not 
reported here, as they are identical to those of Tables [5] and This is because they 
consist of untwisted flavours, the mass tuning of which is independent of K cr . 

There is a rough check which enables us to "predict" the discrepancy in our 
results, induced by an offset 5k ct in the critical point. The standard PCAC depen- 
dence of the squared pseudoscalar mass on the average valence quark masses implies 
that: 

) 2 oc2 am R , w , (A.9) 





(aM«f) cc 2[a^ R , tw + Z m (5am CT )} , (A. 10) 



2 



aM s e d oc am RjW + a/i R)tw + Z m (5am CT ) . (A.ll) 



The above expressions have been obtained by keeping track of the offset in the 
various quark masses (cf. Eqs. (12. 2p . (|A.4h and (|A.7h ) in the tuning procedure, 
through straightforward lowest order Taylor expansions in 5am CT . Since, in the 
absence of offset 5am cr , the quark masses are tuned to satisfy am Rw = a/x Ri t w , we 
easily derive 

gjggV i . Z m (5am CT ) 

aMjg J 2am RjW ' 1 ' ' 

(aMf d \ 2 i Z m (Sam cv ) 

\aMg ) am RjW ' 1 ' ' 
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This means that the offset of the first ratio, due to Sk ci , is "predicted" to be half 
of that of the second. This is roughly confirmed by the data in the case of (3 = 6.1 
with the offset k ct . The relevant results are gathered in Table [TH Clearly, as the 
whole procedure does not take into account higher order discretization effects, our 
expectations are confirmed at a qualitative level. 

We now turn to the 7r/4 case, in which the mass tuning proceeds in a different 
way (cf. ref. [1]). For a fixed bare twisted mass a//Q,tw> the mass degeneracy condition 
n^R,tw = A*R,tw fixes the subtracted mass m qi t w , in terms of Eqs. (|2.5p and (12. 6|) . to 
the value 



™ q ,tw = 7^-a/«o,tw ■{ I + 



ZZ A 



ipa ~ 6m) - ZZkK 



afioMf ■ ( A - 14 ) 



Now this value of am qi t w induces fixing the hopping parameter to say, k or k, 
depending on whether we are working with 

^cr or fc cr • 

When we perform simulations at hopping parameter R (based on the tuning with 
k' ct ), we are not really at subtracted quark mass am qjtw , but rather at am q)tw = 
1/(2k) — 1/(2k ci ). This implies the following offset in the subtracted quark mass: 

arn q ,tw = am q>tw + 5am cr . (A. 16) 

The true untwisted quark masses of our simulation are then 

amR,tw = amR,tw - Z m 8am cr H , (A. 17) 

a#R,tw = a^R,tw H , (A. 18) 

with higher orders omitted. Combining these two expressions, we finally arrive at 
the estimate for the twist angle 

a/^R,tw 

This is again the source of a significant deviation from the target twist angle vr/4; 
the Ward identity results of Table [T^] discussed above corroborate this conclusion. 

In Table [T5l we list the results for the pseudoscalar masses and decay constants. 
They are significantly different to the ones obtained with the new k ct {(5 = 6.1) (cf. 
Table [7|). We also notice that the excellent agreement between ro-Mff and r^M^d m 
Table[7]is lost in Table fl5l Recalling that the two quantities in question, essentially 
being the two sides of a Ward identity, are equal up to discretization effects, we 
interpret this discrepancy as a signal of 0{a 2 ) flavour symmetry violations. The 
comparison of roF S( i to tqF S( i confirms these conclusions. 



25 



The previous analysis of the mass offsets in the tt/4 case suggests two rough 
checks of the observed discrepancies. First, we note that the twist angle "prediction" 
of Eq. (1A.19|1 gives, for the three twisted bare masses used, 

cot (a) ~ 1.50 , (A.20) 
cot (a) ~ 1.54 , (A.21) 
cot (a) ~ 1.64 , (A.22) 

which is in good qualitative agreement with the Ward identity estimates listed in 
Table [15j Second, we compare the pseudoscalar effective masses M|J, computed 
with our hopping parameter k (tuned with n cr ) and listed in Table [7J to the ones 
computed with the hopping parameter R (tuned with k' ct ) of Table [T5l Henceforth, 
the latter quantities are denoted as M^. PC AC suggests that 

aM s e f oc 2aM Rjtw , (A.23) 
aM s e f oc 2aM Rjtw = am\ M + ap,^ tw = 2aM Ritw 

A glance at Table [16] shows that the agreement between the measured quantity 
{aM^/aM^ — l) and the "predicted" value Z m 5am cv /2a^,tw exceeds expecta- 
tions. 



Zj m oam c 



2au,-o 



.(A.24) 
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Appendix B Tables 





(f) 3 - I 


a 
2r 


L 

2r 


K w 


{Htw, «W),tw) 


-^"conf 


6.0 


16 3 x 48 


0.0931 


1.49 


0.1335 


(0.135169,0.03816) 


402 










0.1338 


(0.135178,0.03152) 


398 










0.1340 


(0.135183,0.02708) 


402 










0.1342 


(0.135187,0.02261) 


400 


6.1 


24 3 x 56 


0.0789 


1.89 


0.1343 


(0.1356465, 0.031711) 


100 










0.1345 


(0.1356510, 0.027123) 


100 










0.1347 


(0.1356560, 0.022523) 


122 


6.2 


24 3 x 64 


0.0677 


1.63 


0.1346 


(0.1357800,0.0283240) 


200 










0.1347 


(0.1357825,0.0259850) 


201 










0.1349 


(0.1357866,0.0212897) 


214 


6.3 


24 3 x 72 


0.0587 


1.41 


0.1348 


(0.1358118,0.0246230) 


200 










0.1349 


(0.1358139,0.0222430) 


205 










0.1350 


(0.1358157,0.0198558) 


175 










0.1351 


(0.1358174,0.0174640) 


201 



Table 1: The parameters of the run at twist angle a = n/2. The dataset at f3 = 
is a new run, while all other data are those of ref. [1]. 
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(kf x T 

v a ' a 


a 
2r 


L 

2r 




O/^O.tw 


^conf 


6.0 


24 3 x 48 


0.0931 


2.24 


0.134739 
0.134795 
U.io4ozo 


0.010412 
0.009142 
U.UUooy 1 


200 


6.1 


24 3 x 60 


0.0789 


1.89 


0.135320 
0.135358 
U. lo04Uo 


0.00810 
0.00720 
U.UUolO 


196 


6.2 


32 3 x 72 


0.0677 


2.17 


0.135477 
0.135539 


0.007595 
0.006125 


73 


6.3 


32 3 x 72 


0.0587 


1.88 


0.135509 
0.135546 
0.135584 


0.0076 
0.0067 
0.0058 


76 


6.45 


32 3 x 86 


0.0481 


1.54 


0.135105 
0.135218 
0.135293 


0.01459 
0.01185 
0.01002 


105 



Table 2: The parameters of the run at twist angle a = n/4. The dataset at (3 = 6.1 
is a new run, while all other data are those of ref. [1]. 



& 


K cr 


ref. 


6.0 


0.135196(14) 


[12] 


6.1 


0.135496, 


[3] 




0.135665(11) 


this work 


6.2 


0.135795(13) 


[12] 


6.3 


0.135823 


ZeRo Coll. 


6.4 


0.135720(9) 


[12] 


6.45 


0.135701 


[3] 



Table 3: The critical hopping parameter k ct at several (3- values, obtained from 
several sources. At (3 = 6.1 a second estimate has been computed in the present 
work. When available, errors have been quoted. 
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P a/i Ritw «MR,tw am Ritw cot (a) cot (a) 



Eq. j2SD Eq. (I2~T7D Eq. (12~T1D Eq. (l2~T5l) Eq. (l2~T8l) 



6.0 


0.07292 
0.06016 
0.05176 
0.04321 


0.07458(2) 
0.06129(2) 
0.05245(2) 
0.04360(2) 


0.00427(15) 
0.00255(15) 
0.00141(16) 
0.00026(16) 


0.059 (2) 
0.042 (3) 
0.027 (3) 
0.006 (4) 


0.054 (2) 
0.039 (2) 
0.025 (3) 
0.004 (4) 


6.1 


0.06097 
0.05216 
0.04332 


0.06222(2) 
0.05305(2) 
0.04395(2) 


0.00269(13) 
0.00176(13) 
0.00122(12) 


0.044 (2) 
0.034 (2) 
0.028 (3) 


0.041 (2) 
0.031 (2) 
0.026 (3) 


6.2 


0.05470 
0.05018 
0.04112 


0.05545(1) 
0.05080(1) 
0.04152(1) 


0.00043(09) 
0.00025(09) 
-0.00017(08) 


0.008 (2) 
0.005 (2) 
-0.004 (2) 


0.006 (2) 
0.003 (2) 
-0.006 (2) 


6.3 


0.04769 
0.04308 
0.03846 
0.03382 


0.04833(1) 
0.04364(1) 
0.03890(1) 
0.03418(1) 


0.00109(07) 
0.00100(08) 
0.00072(09) 
0.00054(08) 


0.023 (2) 
0.023 (2) 
0.019 (2) 
0.016 (2) 


0.021 (2) 
0.022 (2) 
0.017 (2) 
0.015 (2) 


6.0 


0.01986 
0.01744 
0.01602 


0.020091(3) 
0.017646(3) 
0.016210(3) 


0.01959(11) 
0.01734(11) 
0.01601(11) 


0.980(6) 
0.989(7) 
0.994(7) 


0.958(5) 
0.968(6) 
0.974(7) 


6.1 


0.01555 
0.01383 
0.01181 


0.015718(2) 
0.013974(2) 
0.011938(2) 


0.01482(10) 
0.01317(10) 
0.01120(11) 


0.948(3) 
0.949(7) 
0.945(9) 


0.933(6) 
0.934(7) 
0.931(9) 


6.2 


0.01465 
0.01182 


0.014780(2) 
0.011922(2) 


0.01540(8) 
0.01261(8) 


1.047(5) 
1.064(7) 


1.034(5) 
1.051(7) 


6.3 


0.01470 
0.01296 
0.01122 


0.014809(2) 
0.013056(2) 
0.011303(2) 


0.01436(9) 
0.01266(9) 
0.01091(10) 


0.973(6) 
0.973(7) 
0.970(9) 


0.963(6) 
0.964(7) 
0.960(8) 


6.45 


0.02823 
0.02294 
0.01940 


0.028417(2) 
0.023085(2) 
0.019523(1) 


0.02733(5) 
0.02203(5) 
0.01852(5) 


0.962(2) 
0.956(2) 
0.950(2) 


0.951(2) 
0.946(2) 
0.941(2) 



Table 4: Renormalized quark masses os/Ur tw an d auiR,tw and twist angles, computed 
from different procedures, indicated by the corresponding equations. Upper table: 
data for the ir/2 case; lower table: data for the 7r/4 case. 
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a/^o.tw 


r M s f 


r Mjf 


r Mf d 


r Mf d 


D.U 


U.UoolO 

0.03152 
0.02708 
0.02261 


z.uoy ) 
1.900(7) 
1.771(7) 
1.620(7) 


z.uyz^o j 
1.907(7) 
1.780(6) 
1.635(6) 


1.865(6) 
1.738(5) 
1.594(5) 


1.865(5) 
1.733(5) 
1.587(5) 


6.1 


0.0317110 
0.0271230 
0.0225230 


2.024(6) 
1.858(7) 
1.693(7) 


2.039(6) 
1.874(7) 
1.715(6) 


2.027(5) 
1.861(7) 
1.703(6) 


2.023(6) 
1.860(6) 
1.700(6) 


6.2 


0.0283240 
0.0259850 
0.0212897 


2.080(6) 
1.981(7) 
1.792(7) 


2.079(6) 
1.980(7) 
1.795(7) 


2.062(6) 
1.964(6) 
1.779(7) 


2.060(5) 
1.959(6) 
1.777(6) 


6.3 


0.0246230 
0.0222430 
0.0198558 
0.0174640 


2.042(9) 
1.953(8) 
1.830(11) 
1.711(10) 


2.050(9) 
1.962(8) 
1.839(10) 
1.722(9) 


2.049(9) 
1.961(8) 
1.839(10) 
1.720(9) 


2.047(9) 
1.958(7) 
1.833(9) 
1.720(8) 



Table 5: The pseudoscalar meson effective masses in the tt/2 case. The physical 
kaon in these units is roME- hys = 1.2544. 
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P 




r F sc 


r F sd 


roF ud 


roF ud 


6.0 


0.03816 
0.03152 
0.02708 
0.02261 


0.4279(48) 
0.4203(40) 
0.4123(39) 
0.4006(33) 


0.4543(52) 
0.4478(49) 
0.4370(43) 
0.4269(41) 


0.4833(52) 
0.4753(51) 
0.4629(43) 
0.4495(42) 


0.4818(57) 
0.4732(57) 
0.4617(49) 
0.4491(47) 




* 


0.3851(58) 


0.4097(70) 


0.4317(67) 


0.4325(73) 


6.1 


0.0317110 
0.0271230 
0.0225230 


0.4653(49) 
0.4417(48) 
0.4306(50) 


0.4795(56) 
0.4559(50) 
0.4437(59) 


0.4939(49) 
0.4731(53) 
0.4549(63) 


0.4936(62) 
0.4716(62) 
0.4536(67) 




* 


0.3919(112) 


0.4009(136) 


0.4121(134) 


0.4093(151) 


6.2 


0.0283240 
0.0259850 
0.0212897 


0.4706(54) 
0.4663(49) 
0.4491(49) 


0.4817(58) 
0.4746(48) 
0.4553(53) 


0.4898(59) 
0.4812(49) 
0.4611(53) 


0.4890(66) 
0.4812(58) 
0.4600(63) 




* 


0.4170(144) 


0.4153(158) 


0.4188(156) 


0.4175(182) 


6.3 


0.0246230 
0.0222430 
0.0198558 
0.0174640 


0.4693(82) 
0.4632(58) 
0.4581(56) 
0.4380(56) 


0.4781(83) 
0.4695(64) 
0.4637(63) 
0.4450(57) 


0.4856(82) 
0.4764(68) 
0.4688(63) 
0.4517(58) 


0.4847(95) 
0.4758(80) 
0.4691(75) 
0.4501(69) 




* 


0.4080(135) 


0.4115(141) 


0.4166(143) 


0.4151(169) 






cont. 


limit 


extrap. 




all beta 




0.430(18) 


0.412(19) 


0.401(19) 


0.396(22) 


w/o = 6.0 




0.437(32) 


0.429(35) 


0.424(35) 


0.425(40) 



Table 6: The pseudoscalar meson decay constants in the ir/2 case. Each decay 
constant value at the physical kaon mass roM^- hys = 1.2544 has been obtained by 
linear extrapolation in the corresponding pseudoscalar mass squared (in r$ units) 
and is reported in the lines indicated by an asterisk. The continuum limit results 
are obtained from linear extrapolations in (a/ro) 2 . 
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r n M e ? 




r n M e ? 
u sd 


U SCI 


6.0 


1.326(4) 
1.253(4) 
1.207(4) 


0.3918(33) 
0.3861(33) 
0.3826(33) 


1.319(4) 
1.244(4) 
1.198(4) 


0.3979(37) 
0.3907(37) 
0.3864(38) 




1.2544 


0.3864(33) 


1 O C A A 

1.2544 


n om '7/' o^7\ 

0.3917(37) 


6.1 


1.235(6) 
1.170(6) 
1.088(6) 


0.3938(44) 
0.3895(43) 
0.3844(42) 


1.231(5) 
1.166(5) 
1.084(6) 


0.4050(54) 
0.4006(54) 
0.3961(54) 




1.2544 


n on / /i -I \ 

0.3951(44) 


1.2544 


0.4U64(53) 


6.2 


1.299(6) 
1.182(6) 


0.4110(63) 
0.4044(63) 


1.295(6) 
1.176(6) 


0.4032(68) 
0.3943(66) 




1.2544 


0.4084(62) 


1.2544 


0.4001(65) 


6.3 


1.338(9) 
1.259(9) 
1.175(10) 


0.4113(66) 
0.4051(63) 
0.3987(59) 


1.337(9) 
1.259(9) 
1.176(9) 


0.4160(81) 
0.4092(78) 
0.4029(76) 




1.2544 


0.4048(61) 


1.2544 


0.4090(75) 


6.45 


2.054(10) 
1.848(11) 
1.702(11) 


0.4776(66) 
0.4569(65) 
0.4433(63) 


2.052(9) 
1.847(9) 
1.701(10) 


0.4871(81) 
0.4673(81) 
0.4544(82) 




1.2544 


0.4089(69) 


1.2544 


0.4216(91) 






cont. 


limit 


extrap. 


all beta 
w/o = 6.0 




0.420(6) 
0.419(9) 




0.424(8) 
0.420(12) 



Table 7: The pseudoscalar meson effective mass and decay constant in the 7r/4 
case. The decay constant values at the physical kaon mass r M* s = 1.2544 have 
been obtained by linear interpolation (extrapolation at (3 = 6.45) in r$M^ d . The 
continuum limit results are obtained from linear extrapolations in (a/ro) 2 . 
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P 


' ™min ~max ' 

. 2r ' 2r _ 


1 lvl sd 


-Tj ALPHA 
K. + 


tjLANL 
K. + 


„ALPHA; w/o c A 
K + 


6.1 


[1.34, 3.08] 


2.039(6) 


1.433(11)(10)(15) 


1.403(11)(14)(18) 


1.339(10)(8)(13) 






1.874(7) 


1.380(13)(10)(16) 


1.351(13)(14)(19) 


1.290(12)(7)(14) 






1.715(6) 


1.332(12)(9)(16) 


1.305(12)(13)(18) 


1.245(11)(7)(13) 






1.2544 


1.219(28)(37) 


1.194(28)(42) 


1.140(25)(31) 


6.1 


[1.18,3.55] 


1.233(6) 


1.239(14)(9)(16) 


1.204(13)(12)(18) 


1.126(12)(6)(14) 




[1.18,3.55] 


1.168(6) 


1.218(15)(9)(17) 


1.183(14)(12)(19) 


1.106(13)(6)(15) 




[1.26,3.47] 


1.086(6) 


1.190(17)(8)(19) 


1.156(16)(12)(20) 


1.080(15)(6)(16) 






1.2544 


1.247(12)(15) 


1.211(12)(16) 


1.134(11)(14) 



Table 8: Results for the pseudoscalar mass and the various ratios H+ at (3 = 6.1. 
The upper (lower) part of the table corresponds to the tt/2 (vr/4) case. The error in 
roM s d is statistical. The three errors of each 1Z + ratio are, in order of appearance: 
(i) due to the statistical fluctuations of the correlations; (ii) due to the errors of 
Z A and Z N \ (hi) the total error from the two previous ones. The results of the 
extrapolations to the physical kaon mass value are shown at the bottom of each 
/3-dataset: the first error of the extrapolated 1Z + is that arising from type-(i) errors 
of the fitted values, while the second is from type-(iii) errors. 






B K (a = ir/2) 


B K (a = 7r/4) 


6.0 


0.911(22) 


0.926(13) 


6.1 


0.824(26) 


0.843(13) 


6.1 


0.812(24) 


0.779(13) 


6.2 


0.828(30) 


0.798(14) 


6.3 


0.789(34) 


0.775(15) 


6.45 




0.789(19) 



Table 9: The RGI bag parameter Bk for all lattice spacings. The second (third) 
column displays data from a = tt/2 (a = rr/4) tmQCD. All results are from ref. [1], 
except those at f3 = 6.1 in boldface, which have been obtained in the present work. 
The errors, which incorporate both statistical and systematic effects, have been 
estimated as described in ref. [1]. 
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" ^.min ~max " 
Xq x 

_ 2r ' 2r _ 


n)M 3 d 


•£>ALPHA 


7£LANL 


^ALPHA; w/o ca 


6.0 


[1.30, d.l/J 


2.092(6) 
1.907(7) 
1.780(6) 
1.635(6) 


2.632(32)(19)(37) 
2.732(43)(19)(47) 
2.809(45)(20)(49) 
2.852(46)(20)(50) 


2.435(30)(25)(39) 
2.528(40)(26)(48) 
2.596(41)(27)(50) 
2.635(42)(28)(51) 


2.293(28)(13)(31) 
2.381(38)(14)(40) 
2.443(40)(14)(42) 
2.480(40) (14) (43) 






1.2544 


3.009(oo)(73) 


2.778(ol)(7o) 


2.ol4(58J(o2) 


6 1 




— . uoy^u j 
1.874(7) 
1.715(6) 


Z.'iJU^/ / V"'-' )\~ J 

2.627(32) (19) (37) 
2.697(40)(19)(44) 


2.574(32)(26)(41) 
2.641(39)(27)(48) 


2.457(30)(14)(33) 
2.521(38)(14)(40) 






1.2544 


3.010(83)(94) 


2.95(8)(10) 


2.817(79)(85) 


6.2 


[1.29,3.05] 


2.079(6) 
1.980(7) 
1.795(7) 


2.414(28)(17)(32) 
2.449(30)(17)(34) 
2.623(39)(19)(43) 


2.400(27)(25)(37) 
2.435(30)(25)(39) 
2.608(39)(27)(47) 


2.312(27)(13)(30) 
2.347(29)(13)(32) 
2.513(38)(14)(40) 






1.2544 


2.92(10)(12) 


2.91(10)(13) 


2.80(10)(11) 


6.3 


[1.23,3.00] 


2.050(9) 
1.962(8) 
1.839(10) 
1.722(9) 


2.352(39)(17)(43) 
2.506(45)(18)(49) 
2.489(56)(18)(59) 
2.661(58)(19)(61) 


2.366(39)(24)(46) 
2.522(45)(26)(52) 
2.504(56)(26)(62) 
2.677(58)(27)(64) 


2.284(38) (13) (40) 
2.433(44)(14)(46) 
2.417(55) (14) (56) 
2.584(56) (15) (58) 






1.2544 


2.95(12)(12) 


2.96(12)(13) 


2.85(11)(12) 



Table 10: Results for the pseudo-scalar mass and the various ratios 1Z~ for twist 
angle ir/2. For an explanation of the errors, see the caption of Table El 
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" ~.min -max " 

_ 2r ' 2r _ 


1 U J "SCl 


•^ALPHA 


T^LANL 


^ALPHA; w/o ca 


6.0 


[1.30,3.17] 
[1.30,3.17] 
[1.30,3.17] 


1.326(4) 
1.253(4) 
1.207(4) 


3.202(40)(23)(46) 
3.306(47) (23) (53) 
3.379(54) (24) (59) 


2.880(37) (31) (48) 
2.971 (42) (32) (53) 
3.035(48) (33) (58) 


2.657(34)(16)(38) 
2.739(40) (16) (43) 
2.798(45)(17)(48) 






1.2544 


3.304(46)(53) 


2.968(41)(53) 


2.738(38)(41) 


6.1 


[1.18,3.55] 
[1.18,3.55] 
[1.26,3.47] 


1.233(6) 
1.168(6) 
1.086(6) 


3.199(67)(23)(71) 
3.296(77)(23)(80) 
3.429(91)(24)(94) 


3.112(65)(32)(73) 
3.204(74) (33) (81) 
3.330(88) (35) (95) 


2.909(62) (17) (64) 
2.995(70) (17) (72) 
3.112(83)(18)(85) 






1.2544 


3.160(62)(65) 


3.077(60)(68) 


2.875(57)(59) 


6.2 


[1.35,3.52] 
[1.35,3.52] 


1.299(6) 
1.182(6) 


2.999(58) (21) (61) 
3.127(72)(22)(75) 


2.976(57)(32)(65) 
3.102(72)(33)(79) 


2.828(55)(16)(57) 
2.945(68) (17) (70) 






1.2544 


3.050(62) (65) 


3.026(61)(68) 


2.875(58)(60) 


6.3 


[1.29,2.94] 
[1.35,2.88] 
[1.47, 2.76] 


1.338(9) 
1.259(9) 
1.175(10) 


2.964(64) (21) (68) 
3.056(73)(22)(76) 
3. 164(85) (22) (88) 


2.983(64)(31)(71) 
3.081(73)(32)(80) 
3.204(85) (33) (91) 


2.848(62) (16) (64) 
2.936(70) (17) (72) 
3.040(82)(17)(84) 






1.2544 


3.065(72)(76) 


3.093(72)(79) 


2.945(69)(71) 


6.45 


[1.30,2.84] 
[1.25,2.88] 
[1.30,2.84] 


2.054(10) 
1.848(11) 
1.702(11) 


2.262(41) (16) (44) 
2.384(53) (17) (56) 
2.484(66) (18) (68) 


2.285(41) (23) (47) 
2.408(53) (24) (58) 
2.508(66) (25) (70) 


2.203(40(12)(42) 
2.322(52)(13)(54) 
2.419(64)(13)(66) 






1.2544 


2.700(93)(98) 


2.725(93)(99) 


2.630(90) (93) 



Table 11: Results for the pseudo-scalar mass and the various ratios TZ- for twist 
angle it J A. For an explanation of the errors, see the caption of Table [H 
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aMo,tw 


cot (a) 


cot (a) 




Ea (12 151) 


Ea (12 181) 


0.0317110 


0.102(2) 


0.103(2) 


0.0271230 


0.141(3) 


0.143(3) 


0.0225230 


0.188(3) 


0.191(3) 


0.00810 


1.435(6) 


1.413(6) 


0.00720 


1.495(7) 


1.473(6) 


0.00615 


1.583(8) 


1.561(8) 



Table 12: Twist angles, computed with an offset critical hopping parameter k' ct ((3 = 
6.1) and from two different procedures, indicated by the corresponding equations. 
The upper (lower) part of the table corresponds to the tt/2 (7r/4) case. 





r Mf 




roMf d 


roF sd 


roF ud 


roF ud 


0.0317110 
0.0271230 
0.0225230 


1.979(6) 
1.812(7) 
1.647(7) 


1.910(6) 
1.739(7) 
1.575(6) 


1.904(6) 
1.735(6) 
1.568(6) 


0.4697(55) 
0.4462(48) 
0.4311(55) 


0.4745(55) 
0.4509(49) 
0.4284(56) 


0.4746(60) 
0.4502(57) 
0.4284(60) 



Table 13: Pseudoscalar effective masses and decay constants, for the tt/2 case, com- 
puted with the offset critical hopping parameter n' cr ([3 = 6.1). 



(3 o^o,tw 


(Mf/M s f) 2 -1 


{MfjM^f-l 


Z m (5am cr .)/am R)W 


6.1 0.0317110 


-0.044 (1) 


-0.109 (2) 


-0.06 


0.0271230 


-0.050 (1) 


-0.124 (3) 


-0.07 


0.0225230 


-0.053 (2) 


-0.134 (4) 


-0.09 



Table 14: Comparison of the flavour breaking effects in the pseudoscalar masses, 
computed with the offset k' cj .((3 = 6.1). In the last column, we also show the results 
of the "prediction" based on Eq. (|A.13p , 
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QMo.tw 


r Q Mf 


r Mf 




roF sd 


0.00810 
0.00720 
0.00615 


1.382(5) 
1.326(5) 
1.257(5) 


1.372(5) 
1.316(5) 
1.247(5) 


0.4001(36) 
0.3944(36) 
0.3862(36) 


0.3325(34) 
0.3196(33) 
0.3025(33) 



Table 15: Pseudoscalar effective masses and decay constants, for the 7r/4 case, com- 
puted with the offset k' ci ((3 = 6.1). 



p 




(Mf d /Mf d ) 2 -1 


Z m (5am cr )/2afj 1 -E i ,w 


6.1 


0.00810 


0.25 (2) 


0.25 




0.00720 


0.28 (2) 


0.28 




0.00615 


0.34 (2) 


0.33 



Table 16: Comparison of the effect of K r cr (f3 = 6.1) on the the pseudoscalar mass, 
measured in the simulation, with the "prediction" based on Eq. (|A.24|) . 
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